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Summary 

When planning an experimental investigation, we are frequently faced with factors that are 
difficult or time consuming to manipulate, thereby making complete randomization impracti- 
cal. A split-plot structure differentiates between the experimental units associated with these 
hard-to-change factors and others that are relatively easy-to-change and provides an efficient 
strategy that integrates the restrictions imposed by the experimental apparatus. Several indus- 
trial and scientific examples are presented to illustrate design considerations encountered in 
the restricted randomization context. In this paper, we propose classes of split-plot response 
designs that provide an intuitive and natural extension from the completely randomized con- 
text. For these designs, the ordinary least squares estimates of the model are equivalent to the 
generalized least squares estimates. This property provides best linear unbiased estimators 
and simplifies model estimation. The design conditions that allow for equivalent estimation 
are presented enabling design construction strategies to transform completely randomized 
Box-Behnken, equiradial, and small composite designs into a split-plot structure. 
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1 Introduction 


Response surface methodology (RSM) is widely used to characterize and optimize products 
and processes. One reason that RSM enjoys success is its recognition of the unique op- 
portunities and restrictions in the industrial context. For example, a sequential strategy is 
employed where small exploratory experiments are used to guide the design of larger more 
comprehensive experiments. Motivated by these opportunities, Box and Wilson [1] adapted 
agricultural methods and developed new experimental strategies tailored to the industrial 
context. In a similar way, split-plot designs provide an efficient statistical method that rec- 
ognizes the common occurrence of hard- or costly-to-change factors resulting in a practical 
restriction on randomization. 

In a split-plot design (SPD), we employ two layers of randomization to reduce the number 
of times that the hard-to-change (HTC) factors are reset. This is a generalization of a 
factorial experiment in which we have two types of experimental units. First, the whole-plot 
factors are randomly assigned to the whole-plot experimental units, which with replication 
generates the whole-plot error variance, of. Then, the subplot factors are randomly assigned 
to experimental units within each whole-plot and generates the subplot error variance, of. 
For a more detailed development of industrial split-plot designs see Montgomery [2] and 
Myers and Montgomery [3]. 

In traditional response surface designs, all factors are assumed to be equivalent in terms of 
ease of manipulation. In this paper, we consider experimental situations where this is not 
the case and explore split-plot designs that accommodate the limitations of the experimental 
apparatus or process. This is similar to the practice of blocking, designed to reduce the effects 
of variability associated with nuisance factors. A block consists of a group of experimental 
runs executed under relatively homogenous conditions. Blocking boundaries, such as days 
or batches, are naturally determined by the restrictions of the experimental process. A 
split-plot design is similar in structure to a blocked design if we view the whole-plots as 
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blocks. However, in a split-plot design we randomly assign factor combinations to the whole- 
plots enabling the estimation and testing of whole-plot effects, in contrast to reducing or 
eliminating them. 

Often there exists a natural whole-plot size determined by the experimental apparatus or 
process, denoted by n. For example, if a temperature chamber can accommodate four subplot 
experimental units, then it would be natural to utilize a design with a whole-plot size of four 
resulting in n min = n max — n — 4. If the hard-to-change factors are truly hard to change 
(in contrast to inconvenient to change), then a minimum number of subplot runs, n m j n , 
should be executed for each setting of the HTC factors. When we evaluate the size of a 
split-plot design, we consider the two types of experimental units and their respective costs. 
For example, if the subplots are significantly less expensive than the whole-plots, then the 
overall cost of the experiment is dominated by the number of whole-plots rather than the 
number of subplot runs. Finally, if n mm = 1, then the completely randomized context is 
more appropriate. 

2 Illustrations of Restricted Randomization 

The following examples illustrate the diversity of industrial applications and motivate the 
number of factors encountered in practice. In Vining, Kowalski, and Montgomery [4] a 
sintering process for the nickel battery plates used in a nickel-hydrogen battery is discussed. 
In this application, the furnace temperature and belt speed are HTC factors. Once these 
factors are set, two composition factors, the density of nickel powder and the amount of 
binder, are relatively easy-to-change (ETC). Once the furnace is stabilized at a specified 
temperature and belt speed, many samples with different composition levels are processed. 
In addition, Vining, Kowalski, and Montgomery [4] analyze an experiment to determine the 
effects of two HTC and two ETC factors on the strength of ceramic pipe. 

Letsinger, Myers, and Lentner [5] discuss a chemical processing experiment to study the 
effect of five factors on a measure of performance. In this experiment temperature and 


3 



pressure are HTC, while the other three factors are relatively ETC. A central composite 
design in five factors was executed in a split-plot structure. In addition, some factors could 
not be set to the specified design points, although their levels were measured and utilized in 
the analysis. This example illustrates a need for a design with performance that is robust 
to set-point deviations. 

Simpson, Kowalski, and Landman [6] consider a wind tunnel experiment involving a race 
car. They study the impact of four factors (front ride height, rear ride height, yaw angle, and 
grille configuration) on the aerodynamic performance. In this study, the front and rear ride 
height factors are hard-to-change because they required considerable time to manipulate. In 
contrast, the yaw angle and grille configuration were easy to manipulate. A second-order 
regression model was developed from the wind tunnel tests that could be used to optimize 
the vehicle performance for a given set of track conditions. 

Micol [7] discusses an experiment to characterize the yawing moment (response) produced by 
the vertical tail configuration of a fighter aircraft as a function of five aerodynamic variables. 
The factors investigated in a wind tunnel experiment are Mach number, Reynolds number, 
angle of attack, angle of sideslip, and rudder deflection. Changing the rudder deflection is a 
manual process requiring a costly and time consuming personnel access into the pressurized 
test section, on the order of 1.5 hours, thereby making it a HTC factor. In addition to 
the time, the facility power requirements to re-stabilize the tunnel conditions after each 
manipulation of this factor represent a significant cost. Alternatively, the other four factors 
can be set remotely within 15 minutes, and therefore are relatively ETC. In this example, 
we see that one setting of the HTC factor is equivalent to 6 settings of the other four factors 
from a time perspective. An experimental strategy that minimizes the number of settings of 
the HTC factor is required. 
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3 Second-Order Split-Plot Designs 


Letsinger, Myers, and Lentner [5] emphasized the need for second-order split-plot designs 
and focused on analysis issues; recommending the use of restricted maximum likelihood for 
model estimation. Draper and John [8] considered modification of the central composite 
design (CCD) (Box and Wilson [1]) to be executed in a split-plot mode to achieve near- 
rotatability resulting in a CUBE or STAR structure. Trinca and Gilmour [9] studied multiple 
layers of split-plotting, called multistratum designs, where the two-level split-plot design is a 
special case of this general class. They employed a computer-intensive search algorithm to 
construct designs that satisfied a minimax criterion of the parameter estimate variance while 
maintaining near-orthogonality between strata. Goos and Vandebroek [10, 11] proposed the 
use of a point exchange algorithm to construct "D-optimal split-plot designs, which minimize 
the volume of the joint confidence region of regression coefficients. Goos [12] provides a 
thorough development of the D-optimal split-plot design approach including implementation 
strategies. 

Vining, Kowalski, and Montgomery [4], hereafter referred to as VKM, proposed a class of 
second-order SPD’s where ordinary least squares (OLS) is equivalent to generalized least 
squares (GLS) for model estimation. In addition, VKM proposed augmentation of these 
designs to provide estimates of pure-error at both the whole-plot and subplot levels. These 
equivalent estimation designs offer numerous benefits. For example, they can be constructed 
such that the equivalence property is independent of the model and variance components. 
Moreover, the parameter estimates are BLUE (best linear unbiased estimators), independent 
of the variance components, and robust to the assumption of normality. With equivalent esti- 
mation designs, model estimation is simplified using OLS, which is available in all statistical 
software packages. 

Parker, Kowalski, and Vining [13] provided a generalized derivation of the equivalence con- 
ditions that leads to the development of systematic construction strategies and the ability 
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to numerically verify that an arbitrary design achieves the equivalence property. They illus- 
trated the construction of split-plot versions of the central composite design. In addition, 
they demonstrated the utility of numerical verification with a class of equivalent estimation 
P-optimal designs. 

In this paper, we review the conditions under which ordinary least squares is an appropriate 
method of estimation for a split-plot design. From these conditions, we illustrate two sys- 
tematic design construction strategies for a split-plot Box-Behnken (Box and Behnken [14]) 
design and discuss cases where a crossed construction is useful. In addition, the ability to 
numerically verify the equivalence property is demonstrated for small composite designs. 

We focus on balanced designs where the whole-plots are of equal size. The experimental 
context determines whether an unbalanced design, with different whole-plot sizes, is a viable 
option. For example, if the subplot experimental units are not significantly less expensive 
than the whole-plots, then the reduction in the total number of subplot runs is desirable. 
However, we believe that in most cases balanced designs are more attractive to the practi- 
tioner, straightforward to execute, and may help to ensure homogeneity among whole-plot 
experimental units. 

4 Equivalence Conditions 

We define the general form of the model as 

y = X/3 + <5 + e , (1) 

where y is the N x 1 vector of responses, X is the N x p model matrix with rank of p, (3 is 
the p x 1 vector of coefficients, S is the N x 1 vector of random whole-plot errors, e is the 
N x 1 vector of random subplot errors, N is the total number of subplot runs, and p is the 
number of terms in the model including the intercept. We assume that S + e has a mean 
of 0 and variance-covariance matrix S = of I + crj J&, where of is the subplot error variance, 
and of is the whole-plot error variance. In the balanced case, Jf, = I m (8) J n , where m is the 
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number of whole-plots, n is the number of subplot runs within each whole-plot, and J„ is a 
n x n matrix of ones. The matrix J& has the form of 

' lnl'„ 0 • ■ ■ O' 

o i n i' n ••• 0 

J b— . . . 

o 0 • • • l n T n _ 

Examining the variance-covariance structure, reveals that observations in different whole- 
plots are independent, while observations within a whole-plot are correlated. 

Parker, Kowalski, and Vining [13] provide a detailed derivation of the equivalence condition, 
which we now summarize. We assume that there is one level of split-plotting and that the 
design supports the model selected. The necessary and sufficient condition for equivalence 
from McElroy [15] (see also Graybill ([16], p. 209) is 

XF = EX, (2) 

where F is a p x p non-singular matrix. We define F to be 

F = of I + ofK, 


where, 

K = (X'X) -1 X'J fe X. (3) 

It can be shown that the condition for equivalence can be expressed as 

XK = J 6 X. (4) 

The right-hand side of this equation sums the columns of the model matrix within each 
whole-plot. These column totals, or unsealed whole-plot means, must equal the model matrix 
multiplied by K to satisfy the condition of equivalence. Note that the equivalence condition 
is independent of the variance components, of and of. 
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For a complete second-order model, let X be partitioned as 



1 

w Dl 

W Ql 

s Dl 

s Ql 

X = 

1 

W Da 

Wq 2 

Sd 2 

Sq 2 


1 

w Dm 

w Qm 

s Dm 



( 5 ) 


where the first column corresponds to the intercept, Wd, denotes the whole-plot main ef- 
fects and two-factor interactions, Wq. contains the whole-plot pure quadratics, Sd ; includes 
the subplot main effects, subplot by subplot two-factor interactions, and the whole-plot by 
subplot interactions, and Sq. denotes the subplot pure quadratic terms of the i th whole-plot. 


Consider one choice for K( pxp ) 


K 


of the form 

" n O' 

0 nlr) 

0 0 
0 0 
0 0 


O' 

0 

0 
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o' m'o ■ 

0 0 
0 Mw , 
0 0 
0 V s 


( 6 ) 


where n is the number of subplot runs per whole-plot, mo is a vector of length k, Id and Iq 
are (w + x (w + w< ' w ^ 1 - 1 ) and w x w identity matrices, respectively, Mw is a w x k 

matrix, Vs is a k x k matrix, where w is the number of whole-plot factors, and k is the 
number of subplot factors. 

By performing the multiplication in Equation (4), we see that, the conditions to achieve 
equivalence are isolated to the relationships between the within whole-plot column sums of 
the intercept and the whole-plot and subplot pure quadratic terms. 

If the design possesses the equivalence property, then (3 OLS = (3 GLS and ordinary least 
squares (OLS) is an appropriate method of estimation. The OLS estimate of the model 
coefficients is Pols = (X'X) _1 X'y. This expression highlights that the parameter estimates 
from an equivalent estimation split-plot design are independent of the variance components. 
By equivalence to GLS, the variance-covariance matrix for the OLS estimate is (X'S _1 X) _1 . 



5 Systematic Design Construction Techniques 


Defining the elements of K based on the parent completely randomized design (CRD) and 
the construction technique leads to systematic approaches for building equivalent estimation 
split-plot designs. In this section, we illustrate the construction method proposed by Vining, 
Kowalski, and Montgomery [4], subsequently referred to as the VKM method, the minimum 
whole-plot (MWP) method proposed by Parker, Kowalski, and Vining [13], and the crossed 
construction method. 

The notation D (w,k) is used to denote a design with w whole-plot factors and k subplot 
factors. Designs in the range of D(w, k ) with w G {1, 2, 3} and k G {1, 2, 3, 4} are considered. 

5.1 VKM Method Box-Behnken 

For the Box-Behnken (BBD), we restrict our discussion to w + k = f < 5, since the BBD 
for / > 5 is based on a partially balanced incomplete block design and does not produce a 
consistent form of K. However, with a suitable assignment of the whole-plot factors to the 
columns of the CRD version of the BBD, the construction methods discussed can be applied 
directly to produce equivalent estimation designs. 

To construct a design of this type, we begin with a completely randomized BBD in / factors. 
The design is rearranged such that the HTC factors remain constant within each whole-plot 
(WP) and the subplot designs are orthogonal for a first-order plus interaction model. The 
center points, consisting of all factors at their zero level in coded units, are placed in a 
separate whole-plot. To maintain balance, subplot runs are replicated as required within 
each whole-plot. Finally, the base design can be augmented by replicating the whole-plot 
containing all center points to provide a pure-error estimate of the whole-plot variance. 

To minimize the number of whole-plots, we use a single whole-plot with the HTC factors 
simultaneously set at it their center level in combination with the pairs of subplot factors 
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forming 2 2 factorial combinations. The size of this whole-plot, for / < 5, is 

Umax = 2 2 ( 2 ) = 2k ( k ~ X )' 

In the other whole-plots, with the HTC factors at ±1, replicates of the base subplot design 
may be required to balance the design. Using this strategy, the whole-plot size can become 
quite large and may not be practical in certain applications. In these cases additional whole- 
plots of smaller size could be used. Nevertheless, we assume that these minimum whole-plot 
designs are permissible. 

We define K for this structure by setting m 0 = 0 and M w = 0. By performing the 
multiplication in Equation (4) the conditions that must be satisfied are 

nW Di = ll'W Di 
nW Qi = ll'W Qi 

and 


ll'Sq, = S Q V s V i. (7) 

For the first set of conditions, since whole-plot factors are held constant within each whole- 
plot, the column sums of the whole-plot model matrix equals n times the level of the whole- 
plot model term. These whole-plot model conditions are satisfied for any design by the 
split-plot structure. 

To define Vg, we consider the column sums of the pure quadratic terms in each whole-plot 
and the condition specified in Equation (7) and find 

V s = 2(k-l)l k l' k . 

As an example, consider a design for D(l,2). Following the VKM construction method we 
begin with a completely randomized BBD in 3 factors, rearrange the design, and add subplot 
center replicates as required to balanced. The resulting design has m — 4 whole-plots of size 
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n — 4 and is shown in Tabic A.l. To specify the K matrix, we find n — 4 and Vg = 2 I 2 I 2 , 
which satisfies the conditions. 

The VKM construction can be applied when / > 5, however the K matrix takes on a slightly 
different form due to the partially balanced nature of the parent design. For the D(l, 3) and 
D(2, 3) cases, a VKM construction BBD has not been identified. For the D(l,3) case, a 
pseudo-crossed approach can be employed where a BBD in the subplot factors is crossed 
with the three levels of the whole-plot factor. The resulting design contains three whole- 
plots with identical subplot designs and one whole-plot containing all centers, see Table A. 2. 
This strategy can also be applied for the D(2,3) case. Note that in the case of D(3, 3) and 
D(3, 4), the subplot pure quadratics are whole-plot level terms. 

The VKM construction method provides a flexible approach to design construction. The 
equivalency property of the design using the VKM method is independent of the model; it 
assumes only that the design supports the model. With the VKM method, a subplot design 
that contains center points, must contain only center points. 

5.2 MWP Method Box-Behnken Designs 

In this section, we propose a class of equivalent estimation designs that require a minimum 
number of whole-plots. To compute the minimum number of whole-plots required for a BBD, 
we consider all combinations of the WP factors at their three-levels (—1, 0, +1). For w = 1, 
rrimin = 3, and for w > 1 and / < 5, 

rrimin = 4 ^ 2 ) + 2w + L 

For example, the minimum number of whole-plots for the D(2,2) case requires, m m i n = 
4(1) + 4+1 = 9. Note that for / > 5, the BBD is based on a partially balanced incomplete 
block design and the number of whole-plots depends on the assignment of the whole-plot 
factors to the columns of the completely randomized design. The development that follows is 
based on f < 5. The construction of designs for / > 5 is not described; however, a summary 
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table of these designs is provided. 


To generate an equivalent SPD, we begin with the completely randomized design and proceed 
in the manner previously discussed by rearranging and ensuring orthogonal subplot designs. 
However, for the MWP method construction, subplot center runs can be added one-by-one 
to each whole-plot maintaining balance. 

For the MWP class of designs, K is defined by setting Vs = 0 in Equation (6). Then, by 
performing the multiplication in Equation (4) we find the condition that must be satisfied, 
in addition to the conditions involving the whole-plot model matrix, is 

inl'nSQ. = lnm'o + Wq.M W V 1 . (8) 

Consider the types of whole-plots based on the level and combinations of the whole-plot 
(WP) and subplot (SP) factors defined by the following three cases; 

1. All WP factors simultaneously at their center level in combination with SP factors at 
their ±1 and center level. 

2. Two WP factors at their ±1 level simultaneously with the SP factors at their center 
level. 

3. A WP factor at ±1 in combination with the other WP factors at their center level and 
one SP factor at its ±1 level. 

Note that Case (2) is not present for w — 1. 

It can be shown, that the column sum of the subplot pure quadratic terms for Case (1) 
is 4 (k — 1) and Wq ; is a matrix of zeroes. Therefore, mo = 4 {k — 1)1*,. For Case (2), 
the column sum equals zero and Wq. = [ l n l n • • ■ 0 ] , where the position of the 

vectors of ones will change depending on which two WP factors are at their ±1 level. For 
Case (3), the column sum equals 2 and Wq. = [ l n 0 ■ ■ ■ 0 ] where the position of 

the vector of ones will change depending on which WP factor is at its ±1 level. We define 
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Mw to be of the form <frl w l' k . For Case (2), we find 


4 (A: — 1) 


For Case (3), we find 


0 = 6 — 4 k. 


Therefore, we summarize the conditions as 


6 — 4 k if w = 1 

— 4( ~ fc ~ 1 ' > =6 — 4 k if w > 1, k > 1 


As a result of this restriction on 0, MWP method BBD’s are not permissible for D(2,3), 
D(2,4), D(3,3), and D(3,4). In addition, other special cases include D(2, 1) and D(3,l), 
where the MWP BBD is equivalent to the VKM construction BBD. 


We note that the relationship to achieve equivalence is not dependent on the whole-plot size, 
n. Therefore, adding subplot center runs to each whole-plot does not change the column 
sums, and thereby maintains the equivalence property. 

Consider a MWP BBD for D(l,2). We find that rn mtn — 3, n — 5, m 0 = 41 2 , and 
M\y = — 210 The design is shown in Table A. 3. Comparing this design to the VKM 
method BBD for D(l,2), we see that there is one less whole-plot and the the whole-plot size 
has increased by one subplot run. 

Now, consider including another HTC factor to the design requiring a D(2,2) design. For 
this design we find, m m i n — 9, n — 5, m 0 = 41 2 , and Mw = — 21 2 1' 2 , which satisfy the 
equivalence conditions. A D(2, 2) design based on a four factor BBD is shown in Table A. 4. 

The MWP method establishes the lower boundary for the number of whole-plots required 
to conduct a Box-Behnken design in split-plot mode. In addition, its structure allows for 
subplot centers to be included within each whole-plot, and the subplot centers can be added 


13 



one-by-one to each whole-plot maintaining the equivalence property. (Note that this feature 
is not available for D(2,l) and D(3,l), which are identical to the VKM method designs 
and therefore do not allow subplot centers to be added one-by-one.) Finally, a whole- 
plot containing all center points for both the whole-plot and subplot factors violates the 
equivalence condition and is not permissible. 

Table B.l provides a summary of the design size for the VKM and MWP method construction 
split-plot BBD’s for all combinations of 1 < w < 3 and 1 < k < 4, where a design has been 
identified. 

5.3 Crossed Designs 

It is well-known that OLS is an appropriate method of estimation of a first-order plus inter- 
action model for a crossed design with orthogonal subplot designs. Bisgaard [17] discusses 
these designs as a Cartesian product and Taguchi [18] refers to them as product arrays. 
From our derivation, we can show that the equivalent estimation property also holds for 
second-order designs. In the model matrix X, we see that Sq ; = Sq,W Referring to K. let 
V s = 0 and M w = 0 and the equivalence condition becomes CI^Sq = l n m' 0 . Therefore, 
m 0 = I^Sq. Therefore, using identical subplot designs results in an equivalent estimation 
second-order design, regardless of the specific design chosen. 

In many cases, a crossed construction results in an unnecessarily large design. However, 
a crossed construction can be an attractive alternative if the parent designs involved have 
relatively few points. In this section we consider crossing combinations of 3-lcvel designs in 
one factor and equiradial designs in two factors. 

We can achieve the minimum design size equiradial design by using a pentagon augmented 
with a center point to avoid singularity of the moment matrix; center points are essential 
to all equiradial designs. However, the pentagon is saturated for a second-order model, and 
therefore it is not recommended unless dictated by cost. Alternatively, a hexagon augmented 
with a center point is an attractive choice, requiring two fewer design points than a spherical 
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central composite design (CCD) and provides one degree of freedom for lack-of-fit. (Note 
that the spherical CCD in two factors is an equiradial octagon design.) The hexagonal design 
features one factor at three levels and the other at five levels, providing a nice compromise 
in the number of required factor levels compared to the CCD. 

First we consider designs at the whole-plot level. For D(l, k), a 3-level design in the whole- 
plot factors can be crossed with a suitable design in the subplot factors. For D(2, k) a 
saturated pentagon plus center point design requires 6 whole-plots, while the hexagon plus 
center point requires 7 whole-plots. It is straightforward to see for subplot designs that a 
3-level design can be applied for D(tc, 1) and a hexagon plus center point for D(tc, 2). 

As an example, consider a design for D(2,2) based on crossing two hexagonal designs at 
the whole-plot and subplot level shown in Table A. 5. This crossed design features two fewer 
whole-plots than the MWP BBD shown in Table A. 4. However, the subplot size has increased 
from 5 to 7. Even with this increase in the number of subplot runs, the reduction in the 
number of whole-plots may be more attractive if the cost of the experiment is dominated 
by setting the HTC factors. There are no replicates at the whole-plot or subplot level 
included in this base design. One approach is to augment the whole-plots with replicated 
subplot centers. In addition, we might consider replicating whole-plot number 7 (whole-plot 
factors simultaneously at their center level) to obtain a pure-error estimate of whole-plot 
error variance. 

6 Numerically Verified Small Composite Designs 

For some classes of completely randomized response surface designs, a consistent analytical 
form of K does not exist. However, configuring these designs into a split-plot version can be 
performed by following a systematic construction method and then verifying if the design 
possess the equivalence property. Using Equation (3) to directly compute K. we can verify 
if the condition from Equation (4) is satisfied. 
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As an example, consider the small composite designs (SCD) proposed by Hartley [19] con- 
structed to reduce the size of the factorial, or cube, portion of a composite design. Draper 
[20] and Draper and Lin [21] proposed using a subset of columns from a Plackett-Burman 
design (PBD) and recommend column choices that provide designs with a repeated run and 
those with high relative V values in a completely randomized context. 

To illustrate the utility of numerical verification using an SCD, consider a split-plot design 
for D(l, 4). We begin with a five-factor cuboidal SCD based on a PBD factorial that includes 
a repeated run, which is deleted, to obtain a saturated design. A center point is added to 
provide one degree of freedom for lack-of-fit, even though it is not essential in a cuboidal 
design. Applying the MWP method and replicating subplot centers as required to achieve 
balance, we obtain the split-plot design shown in Table A. 6. We compute K( 2 ix 2 i) directly 
for a complete second-order model as 


int 
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where the first row indicates the partitioning of the model terms within K. Notice that 
the first whole-plot in Table A. 6 does not contain an orthogonal subplot design resulting in 
non-zero elements associated with the subplot main effects and two-factor interactions in K. 
Applying Equation (4) we find that the conditions for equivalence are satisfied. Inspecting, 
this form of K and considering the equivalence conditions, we can show that subplot centers 
can be added one-by-one within the whole-plots and maintain the equivalence property. 

An alternative design, shown in Table A. 7, was constructed using a complete 12-run PBD 
factorial portion according to the column selection specified by Draper and Lin [21] that 
provides the highest relative V value. This design also has 3 whole-plots of size 9 and 
features orthogonal subplot main effects. While both of these SCD’s reduce the number of 
unique design points, an MWP central composite design based on a resolution V fraction 
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was proposed by Parker, Kowalski, and Vining [13] with 3 whole-plots of size 9, without 
replication (see Table A. 8). From a cost perspective, all of three of these competing designs 
are equivalent, however, their statistical properties are quite different. 

One criterion to compare these three competing designs is their relative D-efficiency, which is 
a measure of their estimation efficiency. Since the designs have the same number of subplot 
runs and the same allocation among three whole-plots, we do not need to apply scaling 
factor. The relative D-effieiency of the minimal point design in Table A. 6 is 49.8% and the 
complete Plackett-Burman factorial design in Table A. 7 is 55.0% compared to the classical 
CCD in Table A. 8. Even if we increase the subplot size to 10, allowing for subplot level 
replication in the CCD, these large improvements in efficiency are maintained. 

Clearly, the CCD is a better choice from an estimation perspective. In addition, the pro- 
jection properties of the CCD are more attractive due to the partial aliasing in the PBD. 
Alternatively, the smaller factorial portion of the SCD provides the ability to estimate sub- 
plot pure-error with fewer subplot runs than the CCD. As in most practical situations, 
compromises among competing criteria in the design selection process must be considered. 

7 Concluding Remarks 

We have emphasized the need to distinguish between HTC and ETC factors and the rel- 
ative cost of their respective experimental units. Integrating apparatus-based restrictions 
into the experimental strategy provides an efficient approach to conduct response surface 
investigations in the presence of restricted randomization. 

Equivalent estimation split-plot designs provide simplified and robust estimation that is inde- 
pendent of the variance components. Simplified because industrial practitioners can perform 
model estimation using OLS software, in contrast to more sophisticated analysis software 
normally required for split-plot designs. Robust in the sense that the design performance can 
be independent of the model, and the parameter estimates are best linear unbiased estirna- 
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tors. Finally, design augmentation providing replication at both the whole-plot and subplot 
level to enable pure-error estimates of the variance components is straightforward. 

Two systematic construction methods to transform a Box-Behnken design into a non-crossed 
split-plot structure have been presented. Cases where a crossed construction is useful have 
been highlighted with three-level and equiradial designs. In addition, a numerical verification 
strategy has been illustrated with small composite designs. A comprehensive catalog con- 
taining equivalent estimation split-plot versions of many classical response surface designs 
is available at the author’s website. The construction methods coupled with the pre-built 
catalog of designs makes equivalent estimation designs available to practitioners. 
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A Equivalent Estimation Split-Plot Designs 


Table A.l: VKM Method Box-Behnken Design for D(l,2). 
(z\ is a HTC factor, x\ and X 2 are ETC factors) 


Whole-Plot 

Zl 

X\ 

%2 

1 

-1 

-1 

0 


-1 

1 

0 


-1 

0 

-1 


-1 

0 

1 

2 

1 

-1 

0 


1 

1 

0 


1 

0 

-1 


1 

0 

1 

3 

0 

-1 

-1 


0 

1 

-1 


0 

-1 

1 


0 

1 

1 

4 

0 

0 

0 


0 

0 

0 


0 

0 

0 


0 

0 

0 
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Table A. 3: MWP Method Box-Bchnkcn Design for D(l,2). 
(z\ is a HTC factor, x\ and x 2 are ETC factors) 


Whole-Plot 

Zl 

Xi 

X 2 

1 

-1 

-1 

0 


-1 

1 

0 


-1 

0 

-1 


-1 

0 

1 


-1 

0 

0 

2 

1 

-1 

0 


1 

1 

0 


1 

0 

-1 


1 

0 

1 


1 

0 

0 

3 

0 

-1 

-1 


0 

-1 

1 


0 

1 

-1 


0 

1 

1 


0 

0 

0 


22 



, X \ and X2 are ETC factors) 


X 2 

Whole-Plot 

Zl 

^2 

X \ 

X 2 

0 

6 

1 

0 

0 

-1 

0 


1 

0 

0 

1 

0 


1 

0 

-1 

0 

0 


1 

0 

1 

0 

0 


1 

0 

0 

0 

0 

7 

0 

-1 

-1 

0 

0 


0 

-1 

1 

0 

0 


0 

-1 

0 

-1 

0 


0 

-1 

0 

1 

0 


0 

-1 

0 

0 

0 

8 

0 

1 

-1 

0 

0 


0 

1 

1 

0 

0 


0 

1 

0 

-1 

0 


0 

1 

0 

1 

0 


0 

1 

0 

0 

0 

9 

0 

0 

-1 

-1 

0 


0 

0 

-1 

1 

0 


0 

0 

1 

-1 

0 


0 

0 

1 

1 

0 


0 

0 

0 

0 

-1 






1 






0 






0 






0 








Table A. 5: Crossed Hexag 
(zi and Z 2 are HTC factors, a 


e-Plot 

A 

^2 

Xi 


1 

1 

0 

1 

0 


1 

0 

0.5 

VT75 


1 

0 

-0.5 

V07T5 


1 

0 

-1 

0 


1 

0 

-0.5 

-V0/T5 


1 

0 

0.5 

-VT75 


1 

0 

0 

0 

2 

0.5 

Vo. 75 

1 

0 


0.5 

V0T5 

0.5 

V07f5 


0.5 

V0T5 

-0.5 

VV75 


0.5 

VT75 

-1 

0 


0.5 

V0T5 

-0.5 

-V07f5 


0.5 

V075 

0.5 

-V07f5 


0.5 

V075 

0 

0 

3 

-0.5 

V0. 75 

1 

0 


-0.5 

V0/f5 

0.5 

V0/f5 


-0.5 

VT75 

-0.5 

VT75 


-0.5 

VT75 

-1 

0 


-0.5 

VT75 

-0.5 

-V07T5 


-0.5 

V0T5 

0.5 

-V0/T5 


-0.5 

VT75 

0 

0 

4 

-1 

0 

1 

0 


-1 

0 

0.5 

V0T5 


-1 

0 

-0.5 

V07f5 


-1 

0 

-1 

0 


-1 

0 

-0.5 

-V0T5 


-1 

0 

0.5 

-V<175 


-1 

0 

0 

0 


and X 2 are ETC factors) 


Whole-Plot 

Zl 


Xi 

x 2 

5 

-0.5 

-V0. 75 

1 

0 


-0.5 

-VT75 

0.5 

VT75 


-0.5 

-V(h75 

-0.5 

V(h75 


-0.5 

-V(h75 

-1 

0 


-0.5 

-V0T5 

-0.5 

-V0T5 


-0.5 

-V(h75 

0.5 

-VT75 


-0.5 

-V(h75 

0 

0 

6 

0.5 

-V0. 75 

1 

0 


0.5 

-V(h75 

0.5 

Vo/f5 


0.5 

-V0T5 

-0.5 

V0/f5 


0.5 

-V0T5 

-1 

0 


0.5 

-V(h75 

-0.5 

-V(h75 


0.5 

-V(h75 

0.5 

-V(h75 


0.5 

-V0T5 

0 

0 

7 

0 

0 

1 

0 


0 

0 

0.5 

V0/f5 


0 

0 

-0.5 

VT75 


0 

0 

-1 

0 


0 

0 

-0.5 

-V(h75 


0 

0 

0.5 

-V0T5 


0 

0 

0 

0 



Table A. 6: MWP Cuboidal Minimal Point Small Composite Design for D(l,4) 
(zi is a HTC factor, x\ , X 2 , x 3 , and X 4 are ETC factors) 


Whole-Plot 

Zl 

Xi 

X 2 

X 3 

X4 

1 

-1 

1 

1 

-1 

1 


-1 

-1 

1 

1 

1 


-1 

1 

-1 

1 

1 


-1 

1 

1 

1 

-1 


-1 

-1 

-1 

-1 

-1 


-1 

0 

0 

0 

0 


-1 

0 

0 

0 

0 


-1 

0 

0 

0 

0 


-1 

0 

0 

0 

0 

2 

1 

1 

-1 

1 

-1 


1 

-1 

1 

1 

-1 


1 

1 

1 

-1 

-1 


1 

1 

-1 

-1 

1 


1 

-1 

-1 

1 

1 


1 

-1 

1 

-1 

1 


1 

0 

0 

0 

0 


1 

0 

0 

0 

0 


1 

0 

0 

0 

0 

3 

0 

-1 

0 

0 

0 


0 

1 

0 

0 

0 


0 

0 

-1 

0 

0 


0 

0 

1 

0 

0 


0 

0 

0 

-1 

0 


0 

0 

0 

1 

0 


0 

0 

0 

0 

-1 


0 

0 

0 

0 

1 


0 

0 

0 

0 

0 
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Table A. 7: MWP Small Composite Design with Complete Plackett-Burman Factorial for 

D(l, 4) 

(zi is a HTC factor, x\ , X 2 , x 3 , and x\ are ETC factors) 


Whole-Plot 

Zl 

x x 

X2 

X 3 

X4 

1 

-1 

-1 

-1 

1 

-1 


-1 

-1 

-1 

-1 

-1 


-1 

-1 

1 

-1 

1 


-1 

1 

-1 

1 

1 


-1 

1 

1 

1 

-1 


-1 

1 

1 

-1 

1 


-1 

0 

0 

0 

0 


-1 

0 

0 

0 

0 


-1 

0 

0 

0 

0 

2 

1 

-1 

-1 

-1 

1 


1 

-1 

1 

1 

1 


1 

-1 

1 

1 

-1 


1 

1 

-1 

1 

1 


1 

1 

-1 

-1 

-1 


1 

1 

1 

-1 

-1 


1 

0 

0 

0 

0 


1 

0 

0 

0 

0 


1 

0 

0 

0 

0 

3 

0 

-1 

0 

0 

0 


0 

1 

0 

0 

0 


0 

0 

-1 

0 

0 


0 

0 

1 

0 

0 


0 

0 

0 

-1 

0 


0 

0 

0 

1 

0 


0 

0 

0 

0 

-1 


0 

0 

0 

0 

1 


0 

0 

0 

0 

0 
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Table A. 8: MWP Cuboidal Central Composite Design for D(l,4) 
(z\ is a HTC factor, x\ , xi , x 3 , and x 4 are ETC factors) 


Whole-Plot 

Zl 

Xi 

X2 

X 3 

X4 

1 

-1 

-1 

-1 

-1 

1 


-1 

1 

-1 

-1 

-1 


-1 

-1 

1 

-1 

-1 


-1 

1 

1 

-1 

1 


-1 

-1 

-1 

1 

-1 


-1 

1 

-1 

1 

1 


-1 

-1 

1 

1 

1 


-1 

1 

1 

1 

-1 


-1 

0 

0 

0 

0 

2 

1 

-1 

-1 

-1 

-1 


1 

1 

-1 

-1 

1 


1 

-1 

1 

-1 

1 


1 

1 

1 

-1 

-1 


1 

-1 

-1 

1 

1 


1 

1 

-1 

1 

-1 


1 

-1 

1 

1 

-1 


1 

1 

1 

1 

1 


1 

0 

0 

0 

0 

3 

0 

-1 

0 

0 

0 


0 

1 

0 

0 

0 


0 

0 

-1 

0 

0 


0 

0 

1 

0 

0 


0 

0 

0 

-1 

0 


0 

0 

0 

1 

0 


0 

0 

0 

0 

-1 


0 

0 

0 

0 

1 


0 

0 

0 

0 

0 
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B Summary of Box-Behnken Designs 


Table B.l: Summary of Box-Behnken Design Size by Construction Method. 

(w is the number of HTC factors, k is the number of ETC factors, m is the number of 
whole-plots, and n is the number of subplot runs per whole-plot) 


Factors VKM MWP 


w 

k 

m 

n 

m 

n 

1 

1 





1 

2 

4 

4 

3 

5 

1 

3 

4 

12 

3 

13 

1 

4 

4 

24 

3 

25 

2 

1 

9 

2 

9 

2 

2 

2 

10 

4 

9 

5 

2 

3 

10 

12 



2 

4 

10 

8 



3 

1 

19 

2 

19 

2 

3 

2 

20 

4 

19 

5 

3 

3 

19 

4 



3 

4 

21 

4 
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